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A  one  dimensional  mathematical  model  is  presented  for  a  high  temperature  lithium-aluminum,  iron 
disulfide  molten  salt  battery  system.  Multi-physics  transport  phenomena  in  the  electrolyte,  charge  bal¬ 
ances  in  the  solid  phases  and  electrolyte,  complex  multi-step  electrochemical  and  chemical  reactions  in 
the  electrodes  are  described  in  this  model.  The  model  includes  the  effects  of  precipitation  salt  on  active 
area  in  the  electrode  and  the  discharge  capacity.  In  addition  the  model  also  takes  the  change  in  volume 
of  the  active  material  into  account  during  the  electrochemical  reactions  by  incorporating  the  change  in 
porosity  of  the  electrode  and  the  change  in  dimension  of  the  electrode.  The  model  results  are  compared 
to  the  existing  model  in  literature  and  available  experimental  data. 

©  2011  Published  by  Elsevier  B.V. 


1.  Introduction 

With  the  oils  crisis  deepens,  the  demand  of  seeking  alterna¬ 
tive  techniques  for  vehicle  propulsion  to  reduce  the  dependence 
on  the  oil  is  growing.  In  addition,  the  requirement  of  the  cheap 
energy  sources  has  given  impetus  to  the  development  of  high- 
energy-density  lithium  batteries  for  various  applications  including 
electric  vehicles,  energy  conversion  and  storage.  To  meet  the  need 
of  the  design  of  an  electric  vehicle  battery,  the  system  must 
have  high  specific  energy  density,  long  life  and  low  cost.  In  addi¬ 
tion,  the  safety,  reliability,  life  and  availability  of  materials  also 
should  be  taken  into  consideration.  As  high  power  sources,  high 
temperature  batteries  are  served  to  provide  power  for  off-peak 
storage  devices  and  electric  vehicle  propulsion  because  of  their 
main  advantages  such  as  long  life,  fast  activation,  rugged  construc¬ 
tion  and  high  reliability  [1-7],  Most  high  temperature  batteries 
use  lithium  alloys  (e.g.  LiAl  or  Li(Si))  as  anodes,  metal-sulfides 
(e.g.  FeS  or  FeS2)  as  cathodes,  and  molten  ionic  conducting  salts 
(e.g.  LiCl-KCl,  LiCl-LiBr-KBr,  LiBr-KBr-LiF,  LiCl-LiBr-LiF)  as  elec¬ 
trolytes.  Among  the  various  cathode  materials,  FeS2  has  its  inherent 
advantages,  such  as  well-predictable  reactions,  cheap  cost  and  non¬ 
toxic  properties  [8],  However,  it  also  has  some  disadvantages  such 
as  moderate  thermal  stability  (about  above  580  °C),  significant  sol¬ 
ubility  in  molten  salts  and  poor  wetting  properties  by  molten  salts 
[8-10],  With  the  help  of  nanotechnology,  the  nanostructured  FeS2 
has  been  prepared  as  a  cathode  material  resulted  in  improved 
electrochemical  performance  [11],  In  addition,  FeS2  electrode  also 
leads  the  fashion  of  the  thermal-sprayed  electrode  [12,13].  All 
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these  developments  give  an  impulse  to  the  LiAl/LiCl-KCl/FeS2  high 
temperature  system  to  have  more  applications  due  to  its  high  open- 
circuit  potential,  high  theoretical  energy  densities  and  inexpensive, 
non-toxic  and  green  compounds  which  help  this  system  to  be  safe. 

In  view  of  the  promising  future  of  the  LiAl/LiCl— KCl/FeS2  high 
temperature  system,  the  investigation  of  the  global  function  of 
the  complex  multi-step  electrochemical,  chemical  reactions  and 
the  multi-physics  phenomena  occurred  in  the  system  is  of  vital 
importance.  Mathematical  modeling  is  an  effective  method  in 
understanding  the  various  complex  processes  involved  in  the 
system.  The  mathematical  model  for  a  LiAl/FeS  battery  system 
has  been  developed  by  Pollard  and  Newman  [3],  Based  on  the 
mathematical  model,  Bernardi  and  Newman  investigated  the  gal- 
vanostatic  discharge  behavior  of  lithium  (alloy),  iron  disulfide  cells 
[1],  However,  these  models  only  predicted  the  occurrence  of  the 
electrolyte  precipitation  salts  based  on  the  experimental  results. 
In  this  work,  the  rate  expression  of  the  precipitation  salts  and  its 
corresponding  criteria  have  been  presented,  and  the  effects  of  the 
precipitation  salts  on  the  porosity  and  discharge  capacity  are  also 
discussed.  Furthermore,  those  models  only  considered  the  change 
in  porosity  of  the  electrode  as  the  resulting  effect  due  to  the  differ¬ 
ences  in  molar  volume  of  reactants  and  various  products  formed 
during  the  electrochemical  reactions.  The  molar  volume  differences 
between  the  reactants  and  products  actually  also  cause  a  change  in 
dimension  of  the  electrode  in  addition  to  the  change  in  porosity.  A 
mathematical  model  accounting  for  the  porosity  change  during  the 
discharge  process  has  been  developed  for  lithium  ion  battery  [14] 
and  lithium/thionyl  chloride  battery  [15]  systems.  In  both  these 
models  the  changes  in  the  active  area  due  to  porosity  change  as 
well  as  the  dimension  change  of  the  electrode  are  incorporated. 
Unlike  the  lithium  ion  battery,  the  volume  change  or  the  expansion 
of  the  electrode  is  considerably  large  in  the  iron  disulfide  battery 
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Nomenclature 

a  specific  area  per  unit  electrode  volume  (cm-1 ) 

do  initial  specific  area  per  unit  electrode  volume 
(cm-1) 

A  cross-sectional  area  of  the  cell  (cm2 ) 

Cp  specific  heat  (J  g_1  K_1 ) 

Deff  effective  diffusion  coefficient  (cm2  s_1 ) 

F  Faraday’s  constant  (96,487  C  equiv.-1 ) 

h  heat  transfer  coefficient  (W  m-2  K_1 ) 

ioj  exchange  current  density  (A  cm-2 ) 

11  superficial  current  density  in  solid  phase  (A  cm-2 ) 

12  superficial  current  density  in  electrolyte  (A  cm-2 ) 
Japp  applied  current  density  (A  cm-2 ) 

Ji  pore  wall  flux  of  species  i  (mol  cm-2  s_1 ) 

kp  precipitation  rate  constant  (cm3  mol-1  s-1 ) 

1°  initial  negative  electrode  thickness  (cm) 

Lp  initial  positive  electrode  thickness  (cm) 

Lp  the  thickness  of  the  positive  electrode  (cm) 

L°  initial  separator  thickness  (cm) 

Lt  total  thickness  of  the  cell  (cm) 

Lk  latent  heat  of  fusion  of  salt  k  (J  mol-1 ) 
m  mass  of  cell  (kg) 

rij  number  of  electrons  transferred  in  electrode  reac¬ 

tionj 

N  number  of  LiAl  particles  per  unit  electrode  volume 

(cm"3) 

rp  radius  of  /J-LiAl  in  LiAl  pellet  (cm) 

ra  radius  of  cr-LiAl  in  LiAl  pellet  (cm) 

R  universal  gas  constant  (8.3145J  mol-1  K_1 ) 

Rg  grid  resistance  ( £2  cm2 ) 

s,j  stoichiometric  coefficient  of  species  i  in  electrode 
reactionj 

S  swelling  coefficient 

t  time  (s) 

t*  transference  number  of  species  i 

T  absolute  temperature  (K) 

Tair  temperature  of  surroundings  (K) 

T0  initial  temperature  (K) 

Uj0  open-circuit  potential  for  reactionj  relative  to  a  ref¬ 
erence  electrode  (V) 

u  local  electrode  displacement  velocity  (cm  s_1 ) 

v*  convective  molar  average  velocity  (cm  s_1 ) 

V  cell  voltage  (V) 

Ve  molar  volume  of  electrolyte  (cm3  mol-1 ) 

Vk  partial  molar  volume  of  salt  k  (cm3  mol-1 ) 

Vj  partial  molar  volume  of  solid  species  i  (cm3  mol-1 ) 

Xj  mole  fraction  of  species  i 

Greek  letters 

oeaj  transfer  coefficient  in  anodic  direction 

acj  transfer  coefficient  in  cathodic  direction 

Ya  activity  coefficient  of  salt  A 

yLl  mean  molar  activity  coefficient  of  lithium  ion 
e  porosity 

£°  initial  porosity  of  the  negative  electrode 

£°  initial  porosity  of  the  separator 

£p  initial  porosity  of  the  positive  electrode 

£j  volume  fraction  of  solid  species  i 

£p  k  volume  fraction  of  precipitation  salt  k 

n'  0i  02  (V) 

r)j  overpotential  of  reactionj  (V) 

iceff  effective  conductivity  of  the  electrolyte  (S  cm-1 ) 

fiA  chemical  potential  of  salt  A 


creff  effective  conductivity  of  the  solid  phase  (S  cm-1 ) 
0i  solid  phase  potential  (V) 

02  electrolyte  potential  (V) 

Subscripts 
app  applied 

A  salt  A,  e.g.,  lithium  chloride 

B  salt  B,  e.g.,  potassium  chloride 

j  reactionj 

k  precipitation  of  salt  A  or  B 

n  negative  electrode 

p  positive  electrode 

ref  reference  electrode 

s  separator 

Superscript 
eff  effective 


system.  Therefore  the  dimension/volume  change  has  to  be  included 
for  accurate  prediction  of  the  performance  of  the  battery. 

In  this  work,  we  have  attempted  to  extend  the  model  presented 
in  Ref.  [1  ]  for  a  lithium-aluminum,  iron  disulfide  high  temperature 
battery  system  by  including  the  expression  of  the  rate  of  precip¬ 
itation  salts  and  the  change  in  dimension  of  the  electrode.  The 
simulations  are  performed  for  battery  performance  at  high  tem¬ 
perature  (450  °C).  The  model  is  validated  by  comparing  with  the 
experimental  data  available  in  the  literature  [16]  and  the  existing 
model  [1], 


2.  Model  development 

The  model  system  consists  of  two  porous  electrodes:  a  posi¬ 
tive  electrode  (FeS2)  and  a  negative  reference  electrode  (LiAl)  with 
a  separator  between  them  and  a  molten  mixture  of  binary  salt 
(LiCl-KCl)  electrolyte  throughout.  Compared  to  Ref.  [1],  we  did 
not  consider  the  reservoir  of  electrolyte,  because  the  emphasis  of 
this  model  is  to  illustrate  the  precipitation  effects  and  the  volume 
change  of  the  electrode  which  all  occurred  in  the  positive  elec¬ 
trode,  neglecting  the  reservoir  of  electrolyte  do  not  alter  the  results 
of  this  model  and  can  be  considered  as  a  valid  assumption.  For 
the  given  schematic  of  the  cell  in  Fig.  1,  during  the  discharge  pro¬ 
cess  lithium  ions  extract  from  the  negative  electrode,  diffuse  and 
migrate  through  the  electrolyte  and  react  with  the  active  mate¬ 
rial  in  the  positive  electrode  to  form  various  intermediate  products 
before  a  complete  discharge.  These  processes  can  be  described  with 
the  help  of  the  porous  electrode  theory  and  transport  equations 
derived  for  a  binary  salt  mixture  by  Pollard  and  Newman  [3,5], 

Physical  processes  for  a  high  temperature  battery  system 
include  the  diffusion,  convection  and  migration  of  electrolyte 
species,  local  reaction  rates,  electrolyte  precipitation  rate,  poten¬ 
tial  of  the  electrolyte  and  solid  phase,  changes  in  the  porosity  and 
electrolyte  composition  due  to  electrochemical  reactions,  energy 
balance  and  local  electrode  displacement  velocities.  The  governing 
equations  for  a  high  temperature  battery  are  summarized  in  the  fol¬ 
lowing  sections  and  the  modifications  made  to  the  original  model 
[3]  are  explained. 

2.1.  Reactions  in  the  cell 

Initially,  the  negative  electrode  consists  of  /8-lithium-aluminum 
particles  which  are  assumed  to  be  spherical  and  non-porous. 
As  the  discharge  process  proceeds,  /8-LiAl  is  oxidized  to  lithium 
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ions  and  aluminum  and  leads  to  a  layer  of  a-Al  forming  at  the 
solid-electrolyte  interface: 

LiAl  -»  Li+  +  Al  +  e-  (Rl) 

The  discharge  reactions  occurred  in  the  positive  electrode  have 
been  investigated  by  many  researchers.  Among  the  various  dis¬ 
charge  mechanisms,  the  most  accepted  mechanism  was  advanced 
by  Argonne  National  Laboratory,  which  follows  four  electron  multi¬ 
stage  reactions  sequence  [7]: 

2FeS2  +  3Li+  +  3e-->  Li3Fe2S4(Z-phase)  (R2a) 


velocity  and  the  transferences  number  are  t*.+  =  xA/2  and  t*+  = 
xBl 2.  The  convection  shown  in  Eq.  (3)  is  considered  negligible  com¬ 
pared  to  the  change  in  molar  volume  of  the  electrolyte  solution, 
A.v*IAVe  <k  1 ,  which  has  been  shown  as  a  valid  assumption  for  this 
system  [4], 

By  neglecting  the  double-layer  charging  and  applying  a  modified 
Faraday’s  law  [3],  the  mass  transfer  rate  of  species  i,  aj,-,  between 
the  electrolyte  and  the  solid  phase  is  given  by: 


■  _  9eP>! 

h  2-j  njF  vk  9t 


(4) 


Li3Fe2S4  +  0.4Li+  +  0.47e~  -*•  1 .58Li2.2Fe0.gS2 
+  0.84Fe0.875S  (R2b) 

Li2+xFei_xS2  +  Fei_xS  o  Li2FeS2  (X-phase)  (R2c) 

Li2FeS2  +  2Li+  +  2e“  2Li2S  +  Fe  (R2d) 

Besides  of  the  above  electrochemical  and  chemical  reactions,  dur¬ 
ing  the  discharge  process,  a  chemical  precipitation  reaction  can  also 
occur  in  the  anode  and  cathode,  respectively,  as  follows,  in  the  melt 
saturated  with  LiCl-KCl  mixtures: 


where  i  represents  the  species  of  Li  or  K,  j  represents  reaction 
(Rl)  for  the  negative  electrode  and  indicates  individual  reaction 
(R2a)-(R2c)  or  (R2d)  for  the  positive  electrode,  k  is  the  precipita¬ 
tion  of  salt  A  (LiCl)  or  B  (KC1),  Sy  is  the  stoichiometric  coefficient 
of  species  i  for  reaction  j  and  n j  is  the  number  of  the  electrons  for 
reaction  j.  The  physical  meaning  of  last  term  on  the  right  side  of  Eq. 
(4)  is  the  removal  rate  of  species  i  from  the  electrolyte  due  to  the 
precipitation  of  salt  k,  and  it  is  given  by  [17]: 

=  fep  x  |§(xfc  -  (xA)kiSat)  (5) 


Li+  +  Cl-^  LiCl(s)  (R3) 

K+  +  Cl-  *=>  KC1  (s)  (R4) 

2.2.  Porous  electrode  model 

The  molten  electrolyte  and  the  solid  phase  are  considered  as  a 
superimposed  continuum  according  to  the  porous  electrode  the¬ 
ory.  The  material  balances  for  the  two  binary  salts  with  a  common 
ion  in  a  macro  homogeneous  medium  are  given  by  the  following 
equations: 


where  kp  is  the  precipitation  rate  constant,  Vk  is  the  molar  volume  of 
salt  k,  {xA  )k  sat  is  the  saturation  composition.  When  the  mole  fraction 
of  salt  k  is  larger  than  the  saturation  composition,  precipitation 
reaction  occurs. 

Because  xB  =  l  -xA,  combining  Eqs.  (1),  (2)  and  (5)  yields  the 
following  equation: 

3exA/Ve  (  xAi2  eDeff  „  \  v^sLi+jQ1j _ 1  dep,Lici 

3t  “  2F  *  %  A)  njF  Vuc^ir 

(6) 


where 

g=t  +  (^+  +  tJ+)| 


i  sDeff  ' 

)  =aJi 

(1) 

=  aJi 

(2) 

(3) 

Ve  is  the  molar  volume  of  the  electrolyte,  xA  and  xB  are  mole  frac¬ 
tions  of  the  salt  A  (LiCl)  and  salt  B  (KC1),  respectively,  a  is  the  specific 
area  per  unit  electrode  volume,  v*  is  the  convective  molar  average 


The  superficial  current  density  in  the  electrolyte,  i2,  is  given  by  the 
following  expression: 

+  (7, 

where  Ke %  is  the  effective  conductivity  of  the  electrolyte,  <p2  is  the 
electrolyte  potential  and  / iA  is  the  chemical  potential  of  LiCl. 

The  current  density  in  the  solid  phase  is  given  by 

U  =  -cref{  V</>,  (8) 

where  oeH  is  the  effective  solid  conductivity  and  <f>\  is  the  solid 
phase  potential. 
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Table  1 

Coefficients  for  the  open-circuit  potential  of  the  FeS2  electrode. 


Reaction  a,-  (V)  bj  (xlO3  VK-1) 

(R2a)  1.4251  0.4785 

(R2b)  1.208771  0.65142 

(R2c)  (in  Eq.  (14))  1.3389  0.0133 

(R2d)  1.43211  -0.147 


A  charge  balance  for  the  current  and  the  electrochemical  reac¬ 
tion  yields  the  following  equation: 

V  ii  =  -V  •  t2  =  -aij  (9) 

where  aij  is  the  local  transfer  current  per  unit  volume  of  the  elec¬ 
trode  due  to  reaction  j  and  can  be  expressed  by  the  polarization 
equation: 


(10) 


where  r)j  is  the  local  surface  overpotential,  pj  =  </>i  -<p2  -  Uj,  and  a  is 
the  specific  area  per  unit  electrode  volume  which  is  assumed  here 
to  be  proportional  to  volume  fraction  of  solid  phase  i,  Sj 

a  =  a0£i  (11) 

where  a0  is  the  initial  specific  area  per  unit  electrode  volume. 

The  exchange  current  density  for  reaction  j,  a0j  ,  can  be  written 
as 


(12) 


where  Yi =  0-5  (assumed  value)  and  xAref  is  equal  to  the  initial  elec¬ 
trolyte  composition. 

The  theoretical  open-circuit  potential  for  reactions  (R2a), 
(R2b)  and  (R2d)  depends  on  the  temperature  and  is  given 
by: 


The  values  for  the  coefficients  a,-  and  bj  are  shown  in  Table  1  [1  ]. 
The  open-circuit  potential  for  reaction  (R2c)  is  given  by: 

_u  (Qc  -  QR2b)  +  T(bc  -  b R2b)  <?R2c2VFeS2 

R2C  R2b  1  -  ((2.2(44  -  2)/(2Xq  -  0.8))  -  3)  [  s°eSJ  J 


where  Xq  depends  on  the  following  linear  variations  of  temperature, 
<7r2c  is  the  number  of  coulombs  of  charge  passed  for  reaction  (R2c) 
per  unit  volume  of  the  positive  electrode: 

x[,  = -9.240  xlO“5T  + 0.91658  (15) 

<}R2c  =  /  aimc  dt  (16) 

Jo 

All  electrode  reaction  potentials  shown  in  Eq.  ( 1 3 )  are  given  relative 
to  a  two-phase  {a- ft)  LiAl  reference  electrode. 

According  to  the  analysis  in  Ref.  [3],  at  the  fully  charged  state  in 
the  negative  electrode,  the  negative  electrode  consists  of  spher¬ 
ical,  nonporous  structured  /3-LiAl  particles.  When  the  discharge 
process  starts,  the  reaction  (R1 )  first  occurs  on  the  outward  region 
of  the  particle  which  results  in  the  establishment  of  a  layer  of  a- 
LiAl  whose  thickness  gradually  grows  at  the  cost  of  fi- LiAl.  The 
Bulter-Volmer  equation  shown  in  Eq.  (10)  is  modified  by  consid¬ 
ering  the  diffusion  overpotential  for  the  mass  transport  of  lithium 


ions  across  the  a-phase  to  the  solid-electrolyte  interface  and  is 
given  by  [5]: 

aij  =  V  •  i2 

= _ exp((aaFIRT)p)  -  exp((-aaFIRT)r, ) _ 

(^/(47tNr^i0  ref(xA/xA  ref)yLi+))  +  ((exp (aaF/RT)p)/(V  •  i2)lim) 

(17) 

where  yLl  is  the  mean  molar  activity  coefficient  of  lithium  ions 
and  (V  •  1*2  )iim  is  the  diffusion-limited  transfer  current  and  can  be 
expressed  by: 


47TiVrI/a(XI  • ) 

(v-i'2)iim=  (i  -x«)((i/rg)-S(ai/r^))  ^ 

where  (x“  )  is  the  saturated  mole  fraction  of  lithium  in  the  a- 
phase,  xg  is  the  average  mole  fraction  of  lithium  in  the  a-phase;  ra 
and  rp  are  the  radius  of  the  a-phase  and  /J-phase  in  the  negative 
electrode,  whose  expressions  are  given  in  Ref.  [5], 


2.3.  Volume  and  porosity  changes  in  solid  phase 


The  volume  change  effects  during  discharge  process  manifest  in 
two  forms:  change  in  porosity  and  change  in  electrode  dimensions. 
This  is  because  when  the  chemical  and  electrochemical  reactions 
occur,  the  products  formed  on  the  surface  of  the  electrode  particles 
occupy  a  part  of  the  pore  volume  and  push  the  electrode  particles 
away  from  each  other. 

The  change  in  porosity  of  the  electrode  during  the  operation 
of  battery  occurs  due  to  the  differences  in  molar  volumes  of  the 
reactants  and  the  products  in  the  chemical  and  electrochemical 
reactions.  The  volume  fraction  of  solid  phase  i  is  given  by: 


3  £j 
3t 


su% 

njFaiJ 


(19) 


where  e,  and  Vt  are  the  volume  fraction  and  molar  volume  of  the 
solid  species  i,  respectively. 

The  appreciable  molar  volume  differences  between  the  reactant 
and  various  products  formed  at  the  cathode  do  not  only  lead  to 
the  change  of  porosity,  but  also  result  in  the  cathode  swelling.  The 
change  in  electrode  dimension  is  governed  by  the  local  electrode 
displacement  velocity,  u,  which  is  given  by  the  following  equation 
[18]: 


<2°> 

where  e  is  the  porosity  of  the  cathode,  S  is  the  swelling  coefficient, 
an  empirical  quantity  and  can  be  estimated  for  the  particular  bat¬ 
tery  system  based  on  the  design  of  the  electrode  and  operating 
conditions  as  demonstrated  in  Ref.  [18],  Eq.  (20)  shows  the  rela¬ 
tionship  between  the  local  electrode  displacement  velocity  and  the 
porosity. 

The  thickness  of  the  positive  electrode  Lp  is  calculated  as  follows: 

T-“'  (2,) 

2.4.  Thermal  effects 

The  temperature  of  the  battery  changes  during  the  discharge 
process  due  to  reversible  heat  of  the  reaction,  heat  losses  to  the 
surrounding  atmosphere,  joule  heating  and  dissipation  of  latent 
heat  because  of  the  precipitation  reaction. 
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The  energy  balance  of  the  cell  is  given  by: 


where  m  is  the  mass  of  the  cell,  Cp  is  the  specific  heat,  A  is  the  cross- 
sectional  area  of  the  cell,  Ujo  is  the  open  circuit  potential  at  450  °C 
for  reaction  j,  h  is  the  heat  transfer  coefficient  and  I*  is  the  latent 
heat  of  fusion  of  salt  k. 

Mole  fraction  of  salt  A  ( xA ),  solid  potential  (0i),  electrolyte 
potential  (02),  volume  fraction  of  precipitation  salt  k  (ep)c),  vol¬ 
ume  fraction  of  solid  species  i  (e,-),  local  electrode  displacement 
velocity  (u),  the  thickness  of  the  cathode  (Lp)  and  temperature  (T) 
can  be  solved  as  the  dependent  variables.  The  governing  equations 
described  by  Eqs.  (5)-(8)  and  (19)-(22)  along  with  the  correspond¬ 
ing  initial  conditions  and  boundary  conditions  described  in  the 
following  section. 


Dimensionless  distance  of  the  cell 


Fig.  2.  Position  dependence  of  mole  fraction  of  Lid  at  different  discharge  times 
{lapp  =  41.6  mA  cm-2 ). 


2.5.  Boundary  and  initial  conditions 

Fig.  1  shows  a  schematic  of  a  high  temperature  cell,  which 
consists  of  four  boundaries  and  three  regions.  At  the  current  collec¬ 
tor/electrode  boundaries,  the  flux  of  the  species  in  the  electrolyte 


At  x  =  0  and  x  =  LT,  VxA  =  0  (23) 

The  solid  phase  potential  is  equal  to  zero  at  the  negative  current 
collector  and  the  solid  phase  current  density  is  set  equal  to  the 
applied  current  density  at  the  positive  current  collector 

At  x  =  0,  0i  =  0 

(24) 

At  X  =  Lj,  i'i  =Iapp 

The  electrolyte  current  density  is  set  equal  to  zero  at  the  current 
collector/electrode  boundaries 

At  x  =  0  and  x  =  LT,  i2  =  0  (25) 

The  local  electrode  displacement  velocity  of  the  electrode  is  set  to 
be  zero  at  the  separator/positive  electrode  boundary 

u\x=l„+ls  =  0  (26) 

There  is  no  electrochemical  reaction  occurred  in  the  separator.  The 
material  balance  of  electrolyte  composition  in  the  separator  is  the 
same  as  that  in  the  cathode  and  anode  except  that  the  porosity  is  a 
constant  in  the  separator.  At  the  interior  boundaries  all  the  fluxes 
are  assumed  to  be  continuous. 

The  initial  values  are  specified  as  follows: 

At  t  =  0,  m  =  n  (negative  electrode),  s  (separator), 
or  p  (positive  elecrode) 

IxA=x° 

Em  =  £m 
U  =  0 

T  =  T0 
Lp  =  I° 


3.  Solution  procedure 

The  numerical  simulations  were  carried  out  in  COMSOL  Multi¬ 
physics  4.1  [19]  by  using  the  finite  element  method  to  discretize 
the  governing  equations. 

3.1.  Discretization  scheme 

Each  model  region  (negative  electrode,  separator  and  positive 
electrode )  is  discretized  along  the  x  coordinate.  The  variables  solved 
for  in  all  three  regions  are:  {xA,  02  and  T};  0 1  is  solved  for  in  the 
positive  and  negative  electrode.  Volume  fraction  of  the  solid  phase 
(gj),  volume  fraction  of  precipitation  salt  k  (£py)  and  the  local  elec¬ 
trode  displacement  velocity  (u)  are  solved  for  only  in  the  positive 
electrode.  To  evaluate  the  dimension  change  of  the  electrode,  the 
thickness  of  the  anode,  separator,  and  cathode  are  converted  to  the 
form  of  dimensionless,  respectively.  Because  the  dimension  change 
is  only  occurred  in  the  cathode  in  this  work,  the  thickness  of  the 
cathode  (Lp),  as  a  variable  is  calculated  based  on  velocity  (u)  at  the 
current  collector/positive  electrode  interface  at  each  time  step. 

3.2.  Parameters 

The  design  parameters,  other  parameters  and  their  values  used 
in  the  simulation  are  shown  in  Table  2.  The  physical  properties 
of  the  electrolyte  and  the  solid  phases  are  listed  in  Appendix  A.  To 
point  it  out,  the  swelling  coefficient  S  is  a  parameter  that  determines 
the  relative  magnitude  of  the  porosity  change  and  the  electrode 
dimension  change.  The  range  of  S  is  from  0  to  1 ,  we  choose  S = 0.5. 
The  number  of  elements  in  the  negative  electrode  and  separator 
is  30;  the  number  of  elements  in  the  positive  electrode  is  96  with 
1.3  element  growth  rate.  The  relative  tolerance  is  le-6  and  the 
absolute  tolerance  is  le— 4. 

4.  Results  and  discussion 

4.1.  Electrolyte  composition 

Fig.  2  shows  the  electrolyte  composition  profiles  across  the  cell 
sandwich  for  several  times  with  a  constant  discharge  current  den¬ 
sity  of  41.6  mA cm-2.  A  distinct  electrolyte  composition  profile  is 
established  on  the  base  of  dimensionless  thickness.  Before  the  dis¬ 
charge  process  starts,  the  mole  fraction  of  salt  LiCl  remains  as 
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Table  2 

Parameters  used  in  the  simulation. 


Parameter  Cathode 

Anode 

Separator 

Reference 

Electrode  parameters 

L°m  (cm)  0.2146 

0.4048 

0.16 

HI 

£°  0.75 

0.39 

0.75 

[1] 

Parameter 

Cathode 

Anode 

Reference 

Kinetics  parameters 

1 

0.5 

[3] 

"j 

2.0 

1.0 

[3] 

Vj 

-2.0 

1.0 

[3] 

Scrj 

0.0 

0.0 

[3] 

0 KaU 

0.68 

0.68 

[1] 

i0  (A cm-2) 

2.8 

[3] 

(a0i0)j  (Acm-3) 

2.0e3 

[3] 

D„(cm2s-1) 

4.0e-10 

[3] 

*LI 

0.05 

[3] 

yU+ 

0.5 

[3] 

Other  parameters 

<7FeS2  ((£2cm_1)  ) 

100 

[i] 

CtLyFejS^Z-phase)  ((£2  Cm-1)  ) 

0.2 

[i] 

^Li2FeS4 (X-phase)  ((£2 Cm-1)  ') 

5.1 

[i] 

oFe  (£2  cm-’r' 

8.3  e4 

[ii 

<fu2s  (£2 cm-')-1 

0 

[i] 

VFes2  (cm3  moC1 ) 

23.93 

[ii 

Vz-phase  (cm3  mol-1) 

99.17 

[ii 

Vx-phase  (cm3  mol-1) 

49.82 

[ii 

VFe  (cm3  mol-1 ) 

18.49 

[ii 

Vu2s  (cm3  mol-1 ) 

28.062 

[ii 

SFeS2,3o 

-4/3 

[ii 

Sz-phase,3a 

2/3 

[ii 

-2 

[ii 

4 

[ii 

[ii 

SLi2S,3d 

2 

[ii 

1 

iij 

VKCi  (cm3  mol-1 ) 

20.5 

[3] 

t-xci  0  mol-1 ) 

25,500 

[5] 

N(lcm-3) 

0.5e7 

[5] 

mCp/A(  kJkg-'K-1) 

F(Cmol_1) 

R  Q  mol-1  K-1) 

1.89 

96,487 

8.3143 

[5] 

h  (W  irr2  K~' ) 

TairW 

8.25e-2 

298.15 

[3] 

Rg  (£2  cm2) 

1.55 

[3] 

S 

0.5 

Assumed  value 

kp  (cm3  mol-1  s-1 ) 

0.001 

Assumed  value 

a  constant  that  is  the  initial  condition  0.58.  When  the  discharge 
begins,  lithium  ions  are  introduced  into  the  electrolyte  from  the 
negative  electrode  and  transport  across  the  separator  to  the  pos¬ 
itive  electrode,  where  they  can  react  with  FeS2  and  finally  form 
Fe  and  Li2S.  Fig.  2  shows  that  as  the  discharge  proceeding,  the 
mole  fraction  of  lithium  ions  increases  in  the  anode  due  to  reaction 
(Rl).  All  the  discharge  reactions  occurred  in  the  cathode  consume 
lithium  ions,  which  deeply  decrease  the  amount  of  lithium  ions. 
At  the  end  of  the  discharge,  the  mole  fraction  of  lithium  ions  is 
almost  down  to  zero.  A  constant  composition  gradient  is  observed 
across  the  separator  due  to  the  constant  flux  of  lithium  ions  in  the 
separator. 

The  details  of  the  composition  variations  can  be  attributed  to 
the  combined  effects  of  diffusion,  migration,  electrochemical  reac¬ 
tions  and  phase  transitions.  In  the  electrodes,  a  small  change  of 
the  local  porosity  may  intensify  the  variations  of  the  electrolyte 
compositions. 

4.2.  Volume  fractions  of  solid  phases  distribution 

According  to  the  discharge  mechanism  proposed  by  ANL,  there 
are  four  intermediate  phases,  Li3Fe2S4  (Z-phase),  Li2+xFei_xS2, 


U,  Utilization  of  FeS2elctrode  (%) 


Fig.  3.  Volume  fractions  of  solid  phases  on  the  interface  between  the  cathode  and  the 
current  collector  as  a  function  of  utilization  of  FeS2  electrode  ( lapp  =41.6  mA  cm-2). 


Fe1_x/SandLi2FeS2  (X-phase).FortheLi2+xFei_xS2  and  Fei_x/ Scorn- 
pounds,  we  do  not  have  molar  volume  data  because  they  are  the 
functions  of  the  composition  variable  x.  We  simplify  the  phase  tran¬ 
sitions  as  Ref.  [1]  did.  Reactions  (R2b)  and  (R2c)  can  be  combined 
into  one  reaction 

Li3Fe2S4  (Z-phase)  +  Li+  +  e"  ^  2Li2FeS2  (X-phase)  (R5) 

Based  on  this  approximation,  the  intermediate  phases  are  Li3Fe2S4 
(Z-phase)  and  Li2FeS2  (X-phase).  We  calculate  the  volume  fractions 
of  the  solid  phases  in  the  positive  electrode  according  to  three  reac¬ 
tions  (R2a),  (R5)  and  (R2d).  This  simplification  is  accurate  because 
reactions  (R2b)  and  (R2c)  are  the  intermediate  reactions  and  can 
reach  the  completion  at  the  end  of  the  discharge  process.  Fig.  3 
shows  the  volume  fractions  of  the  solid  phases  in  the  positive  elec¬ 
trode.  The  utilization  is  defined  by  the  following  expression  which 
can  be  considered  as  a  dimensionless  time  during  the  discharge 
process: 

%  utilization  of  FeS2  =  ^^apptVres2  (28) 

LPB  FeS24f 

where  VVeSj is  the  molar  volume  of  FeS2,  L®  is  the  initial  thickness 
of  the  positive  electrode,  e°eS  is  the  initial  porosity  of  the  positive 
electrode  and  Iapp  is  the  applied  current  density. 

The  species  considered  to  exist  during  the  discharge  process  in 
the  positive  electrode  include  the  initial  active  material:  FeS2,  the 
intermediate  products:  Z-phase  and  X-phase  and  the  final  prod¬ 
ucts:  Li2S  and  Fe.  Fig.  3  presents  the  volume  fractions  of  solid  phases 
on  the  interface  between  the  cathode  and  the  current  collector  as 
a  function  of  utilization  of  FeS2  electrode.  The  simulation  results 
are  under  the  condition  of  68%  (mole  ratio)  lithium  chloride  in  the 
electrolyte.  The  volume  fraction  of  the  FeS2  decreases  with  the  uti¬ 
lization  of  FeS2  during  the  discharge  process.  The  volume  fraction 
of  Z-phase  first  increases  at  about  9%  utilization  and  then  decreases. 
The  reduction  of  the  volume  fraction  of  Z-phase  indicates  the  rate 
of  reaction  (R2a)  is  slower  than  the  rate  of  reaction  (R2b).  When  the 
volume  fraction  of  Z-phase  falls  to  a  small  and  stable  value,  it  means 
that  the  reaction  (R2b)  approaches  to  the  completion.  The  volume 
fraction  of  X-phase  has  the  similar  trend  as  the  volume  fraction  of 
Z-phase,  but  left  behind,  which  confirms  that  the  sequence  of  the 
reactions  occurred  in  the  positive  electrode  is  correct.  For  the  final 
products  Fe  and  Li2S,  their  volume  fractions  always  increase  as  the 
discharge  process  proceeding.  Compared  to  the  volume  fraction  of 
Fe,  the  amount  of  the  volume  fraction  of  Li2S  is  much  larger.  This  is 
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U,  Utilization  of  FeS2 electrode  (%) 


U,  Utilization  of  the  FeS2electrode  (%) 


Fig.  4.  Comparison  of  simulated  discharge  profile  with  experimental  data  and  model  Fig.  5.  Model  results  of  the  cell  voltage  at  different  current  densities, 

results  in  literature  ( lapp  =  50  mA  cm-2 ). 


firstly  because  the  producing  amount  of  Li2S  is  two  times  of  Fe  due 
to  the  stoichiometric  coefficient  of  reaction  (R2d),  secondly  because 
the  molar  volume  of  Li2S  is  much  larger  than  the  mole  volume  of 
Fe. 

4.3.  Discharge  behavior 

Fig.  4  shows  the  cell  voltage  profile  for  galvanostatic  discharge, 
as  a  function  of  the  utilization  of  the  FeS2  electrode.  As  shown 
in  Fig.  4,  it  is  apparent  to  distinguish  the  cell  voltage  profile  con¬ 
sists  of  four  regions.  The  four  regions  correspond  to  the  four-stage 
discharge  reactions  sequence.  Fig.  4  also  shows  the  comparison 
of  the  discharge  curves  between  the  present  model  results,  the 
model  results  from  Ref.  [1]  and  the  experimental  data  from  Ref. 

[14],  The  value  of  S= 0.5  was  used  to  simulate  the  discharge  profile 
for  a  lithium  rich  electrolyte  (68%  Lid)  at  450  °C  with  a  discharge 
current  density  of  50  mA cm-2.  The  same  simulation  parameters 
were  applied  to  the  experimental  data  and  Bernardi’s  model.  The 
discharge  profile  obtained  from  the  model  presented  in  Ref.  [1] 
did  not  consider  dimension  change  in  the  positive  electrode.  The 
model  without  considering  the  volume  expansion  under-predicts 
the  capacity  of  the  battery  severely  (20%).  It  is  therefore  demon¬ 
strated  that  the  inclusion  of  volume  expansion  predicting  the 
discharge  profile  is  more  accurately.  This  is  because  that  in  the 
experimental  battery  system,  changes  in  volume  of  the  reactants 
lead  to  plugging  of  the  pores  as  well  as  expanding/contracting  the 
electrode.  The  relationship  between  the  expansion  velocity  and  the 
porosity  is  established  in  the  present  study.  This  has  resulted  in  g 

an  increase  in  the  utilization  of  the  electrode  and  thus  agreed  the  a> 

capacity  of  the  battery  better  with  the  experimental  data  compared  -3 

to  the  model  in  Ref.  [1  ].  fe 

Fig.  5  illustrates  the  effect  of  increasing  discharge  current  den-  £ 
sity  on  the  discharge  behavior.  It  is  apparently  observed  that  the  h 
active  material  utilization  is  limited  at  higher  discharge  rates.  At 
the  100  mA  cm-2  discharge  current  density,  the  cell  potential  drops 
sharply  when  about  35%  of  the  cathode  material  is  utilized,  this  is 
because  the  overpotential  induced  by  the  higher  discharge  current 
density  increases,  in  turn,  the  domination  of  the  reactions  (R2b)  and 
(R2c)  are  initiated  sooner  than  for  the  case  with  a  lower  current 
density  (e.g.,  50  mA cm-2).  For  an  even  higher  discharge  current 
density,  for  instant,  200  mA cm-2,  we  even  could  not  observe  the 
four-region  discharge  behavior  because  the  discharge  current  den¬ 
sity  is  so  high  that  hinder  the  utilization  of  the  active  material.  In 


this  case,  before  the  reaction  (R2b)  predominated,  the  cell  voltage 
has  dropped  to  the  cell  cutoff  voltage.  From  Fig.  5,  we  found  that 
the  higher  the  current  density  is,  the  faster  the  voltage  cell  drops 
to  reach  the  cell  cutoff  voltage  1 V. 

4.4.  Temperature  behavior 

As  we  discussed  above,  the  temperature  coefficient  fa;  (shown 
in  Table  1 )  of  individual  reaction  in  cathode  can  give  us  an  indi¬ 
cation  of  the  temperature  behavior  of  a  cell  during  the  discharge 
process.  When  the  temperature  coefficient  fa,-  is  positive,  it  means 
the  reaction  is  endothermic;  consequently,  the  temperature 
would  decrease.  In  contrast,  when  fa,  is  negative,  the  counter¬ 
part  temperature  would  increase.  Fig.  6  shows  the  temperature 
change  on  the  interface  between  the  cathode  and  the  current 
collector  as  a  function  of  discharge  time.  It  is  observed  that  the 
temperature  decreases  sharply  before  the  utilization  reaching  42%, 
this  is  because  before  the  utilization  approaching  42%,  reaction 
(R2a)  with  a  positive  temperature  coefficient  (0.475  x  10-3  VK-1) 
predominates  in  the  discharge  reactions  sequence  and  the  rest 
three  reactions  play  slight  roles  in  the  change  of  cell  temperature. 


Fig.  6.  The  temperature  on  the  interface  between  the  cathode  and  the  current  col¬ 
lector  as  a  function  of  discharge  time. 
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U, Utilization  of  the  FeS2  electrode  (%) 


Fig.  7.  Volume  fraction  distribution  of  precipitation  KC1  in  the  positive  electrode  for 
various  discharge  times  with  kp  =0.001  cm3  mol-1  s_1. 

When  the  utilization  proceeds,  reaction  (R2a)  has  reached  its 
completion,  instead  reactions  (R2b)-(R2d)  contribute  a  lot  to 
the  temperature  change.  In  this  region,  the  temperature  declines 
mildly  firstly  because  of  the  negative  temperature  coefficient  of 
reaction  (R2d)  (-0.147  VK_1),  secondly  because  the  precipitation 
of  KCl  occurs  and  the  polarization  heat  becomes  increasingly 
important.  At  the  end  of  the  discharge  process,  one  reason  why  the 
temperature  increases  a  little  rather  than  continuously  decreasing 
is  that  reaction  (R2d)  becomes  the  dominating  reaction  and  it  is 
exothermic.  Another  reason  responding  for  the  phenomena  is  still 
the  polarization  heat  due  to  the  KCl  precipitation. 

4.5.  Effect  of  precipitation  of  salt  KCl  on  the  cell  behavior 

In  the  anode,  lithium  ions  produced  in  reaction  (R3)  may  react 
with  chloride  ions  in  the  electrolyte  to  make  lithium  chloride  crys¬ 
tals  that  precipitate  out  of  the  electrolyte  to  form  solid  lithium 
chloride,  meanwhile,  in  the  cathode  where  lithium  ions  are  con¬ 
sumed  causing  potassium  ions  exceeded.  The  excess  amount  of 
potassium  ions  may  also  react  with  chloride  ions  to  form  solid 
potassium  chloride  according  to  reaction  (R4).  However,  based  on 
the  analysis  in  Ref.  [3],  the  precipitation  of  Lid  is  avoided,  so  the 
precipitation  of  KCl  is  the  emphasis  in  this  work.  When  the  reac¬ 
tions  occurred  in  the  positive  electrode  consume  lithium  ions,  the 
mole  fraction  of  salt  LiCl  decreases  below  the  saturation  value  of 
(xA)sat,Kco  which  results  in  the  onset  of  KCl  precipitation.  Fig.  7 
shows  the  distribution  of  the  volume  fraction  of  solid  KCl  in  the 
positive  electrode  at  different  discharge  times  when  the  precipita¬ 
tion  rate  constant,  kp,  is  0.001  cm3  mol-1  s-1.  During  discharge,  at 
the  separator-positive  electrode  interface,  the  volume  fraction  of 
solid  KCl  is  zero.  With  shifting  inward  from  the  separator-positive 
electrode  interface,  the  volume  fraction  of  solid  KCl  occurs  at  about 
25%  depth  of  the  positive  electrode.  With  discharge  time  increases, 
the  volume  fraction  of  solid  KCl  increases. 

The  precipitation  of  KCl  salt  on  the  electrode  partially  blocks  the 
interface  between  the  electrode  and  electrolyte,  in  turn  reduces 
the  electrochemical  active  surface,  consequently,  and  causes  the 
change  in  porosity.  In  addition,  the  variation  of  porosity  in  the 
electrode  is  due  to  the  electrochemical  reactions  and  the  parasitic 
side  reaction.  In  the  positive  electrode  the  side  reaction  resulting 
in  the  occurrence  of  KCl  precipitation  also  leads  to  the  change  in 
porosity.  Fig.  8  presents  the  precipitation  effect  on  the  change  in 
the  porosity  at  the  interface  between  the  cathode  and  the  current 
collector  as  a  function  of  utilization  of FeS2  electrode.  The  porosity 


Fig.  8.  Change  in  porosity  on  the  interface  between  the  cathode  and  the  current 
collector  as  a  function  of  utilization  of  FeS2  for  different  KCl  precipitation  rate  con- 


U,  Utilization  of  the  FeS2electrode  (%) 

Fig.  9.  Voltage  change  with  utilization  of  the  FeS2  electrode  at  a  discharge  current 
density  of  50  mA  cm-2  for  different  KCl  precipitation  rate  constants. 


at  the  interface  between  the  cathode  and  the  current  collector 
decreases  with  the  utilization  of  the  FeS2  electrode  increases,  and 
the  decrease  along  with  the  utilization  of  the  FeS2  electrode  in 
the  porosity  with  a  large  kp  is  sharper  than  that  with  a  small  kp. 
As  shown  in  Fig.  8,  the  porosity  does  not  change  much  at  the  end 
of  the  discharge  process  (about  90%  utilization  of  FeS2  electrode), 
when  kp  =  0.001  cm3  mol-1  s-1.  However,  when  kp  increases  to 
0.1  cm3  mol-1  s-1 ,  the  porosity  decreases  sharply  at  50%  utilization 
of  the  FeS2  electrode,  which  may  limit  the  discharge  capacity. 

To  elucidate  the  effect  of  kp  on  battery  performance,  the  change 
in  cell  voltage  along  with  the  utilization  of  the  FeS2  electrode  for 
different  values  of  kp  is  presented  in  Fig.  9.  The  discharge  capacity 
increases  with  decreasing  kp.  From  the  discharge  curve,  we  can  dis¬ 
tinguish  the  four  regions  clearly,  which  supports  the  four-region 
behavior  of  the  FeS2  electrode.  However,  with  larger  value  of  kp 
(e.g.  0.1  cm3  mol-1  s— 1 ),  the  porous  electrode  becomes  plugged  as 
shown  in  Fig.  8,  so  that  the  discharge  process  of  the  FeS2  electrode 
is  incomplete.  In  order  to  avoid  this  situation,  the  precipitation  of 
KCl  should  be  decreased  by  increasing  the  temperature,  the  LiCl 
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U,  Utilization  of  the  FeS2  electrode  (%) 


Fig.  10.  Expansion  of  cathode  as  a  function  of  cathode  utilization  for  different  values 
of  swelling  coefficient  S. 


concentration  in  the  electrolyte,  the  porosity  of  the  electrode  and 
so  on. 

4.6.  Dimension  change  in  the  positive  electrode 

The  change  in  dimension  of  the  positive  electrode  during  the 
discharge  process  is  shown  in  Fig.  10  for  various  values  of  swelling 
coefficients.  It  can  be  observed  from  the  profile  that  the  amount 
of  volume  expansion  increased  with  an  increase  in  the  value  of 
S.  When  S  =  0,  where  all  of  the  volume  change  results  in  changing 
in  the  porosity,  the  dimension  of  the  electrode  does  not  change, 
which  is  the  situation  Bernardi’s  model  considered.  After  extend¬ 
ing  Bernardi’s  model  by  including  the  local  electrode  displacement 
velocity,  the  trend  of  the  dimension  change  in  the  cathode  is  seen 
clearly:  As  S  increases,  the  thickness  expansion  of  the  electrode 
keeps  increasing  with  the  utilization  of  FeS2  electrode. 


5.  Conclusions 

A  mathematical  model  for  a  LiAl/LiCl-KCl/FeS2  cell  has  been 
developed  which  can  be  used  to  predict  the  variation  of  the  con¬ 
centration  of  the  ionic  species  in  the  melt  and  the  volume  fractions 
profiles  of  individual  solid  species  as  a  function  of  time  and  space. 
Discharge  behavior,  the  high  discharge  current  density  effect  and 
cell  temperature  behavior  are  also  investigated.  Moreover,  our 
model  predicts  that  KC1  precipitation  can  limit  the  utilization  of 
the  FeS2  electrode  and  discharge  capacities.  To  enhance  the  capac¬ 
ity,  the  precipitation  of  KC1  should  be  decreased  by  increasing 
the  cell  temperature,  the  LiCl  concentration  in  the  electrolyte, 
and  the  porosity  of  the  electrode.  Furthermore,  our  model  is  also 
developed  to  account  for  the  volume  change  in  high  tempera¬ 
ture  battery  system.  The  volume  change  in  the  cathode  during 
discharge  is  attributed  to  the  change  in  porosity  and  the  change 
in  the  dimension  of  the  electrode.  Based  on  the  material  balance, 
the  relationship  between  the  porosity  and  the  local  electrode  dis¬ 
placement  velocity  is  presented.  It  was  found  that  inclusion  of  the 
change  in  dimension  of  the  electrode  is  vital  for  accurate  prediction 
of  the  discharge  curve.  The  model  agrees  well  with  the  experi¬ 
mental  data  and  could  be  used  to  optimize  the  design  parameters 
for  other  operating  conditions  of  lithium-aluminum,  iron  disulfide 
high  temperature  battery  system. 
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Appendix  A. 


The  physical  properties  of  the  electrolyte  are  given  below  based 
on  Refs.  [3,6], 

The  effective  conductivity  of  the  electrolyte  is  given  by  the  fol¬ 
lowing  expression: 

Keff  (fi  cm"* 1 )  =  £15(18.9803x2 3 4 5 6 7 8 9  -  5.0170x„ 

+  9.0903)  exp  (~1475'76)  (A.1) 

The  chemical  potential  gradient  of  salt  A  is  given  by: 

WptA  0mor1cm-1)  =  RT(l+^^)vinx/,  (A.2) 

where  the  activity  coefficient  is  given  by  the  following  expression: 


In  ya  = 


-1555.11 

T 


tanh 


f  1023.84(1  -xAf\ 

V  1555.11X4  J 


(A.3) 


The  effective  diffusion  coefficient  of  lithium  ions  in  the  electrolyte 
is  given  by: 

DeJf  (cm2  s_1 )  =  2.442  x  10-4  exp  ^  14^5.76  j  ^  +  o.4)e°-5 

(A.4) 


Applying  Maxwell’s  model,  we  assume  that  the  solid  phase  consists 
of  a  continuous  phase  and  several  discrete  phases.  In  our  system, 
Fe  and  voids  are  assumed  to  be  the  discrete  phases  and  other  elec¬ 
trode  phases  are  lumped  together  and  called  the  continuous  phase. 
Consequently,  the  solid  effective  conductivity  is  given  by: 


(g*/Qc)-l 

2  ({otf/fft)+y' 


( 


[crefi /a c)  -  1  Y„ 
2{{aeS fac))-¥  1 )  Fe 


(A.5) 


where  ac  is  the  conductivity  of  the  continuous  phase 


 £Z-phase°Z-pha 


ise  +  £X-phase°X-phase  +  eli2S  +  °Li2S 

£z-phase  +  «X-phase  +  £Li2S  +  eFeS2 


+  £FeS2<7FeS2 

(A.6) 


The  molar  volume  of  the  electrolyte  is  defined  as  the  following 
expression: 


(74.555  -  32.161  x  xA)  x  10“16 

e  (2.16233194  -  0.00073684)  -  0.15  x  (x„  -  0.00147368  x  T  -  0.36253409) 

(A.7) 

The  equilibrium  composition  that  characterized  the  onset  of  pre¬ 
cipitation  of  KC1  is  approximated  by: 

(*a)KCi,sat  =  0.14615  +  1.9117  x  10_3T  -  1.948  x  10_6T2  (A.8) 
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